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ABSTRACT 

We cross-correlate foreground cleaned Planck Nominal Cosmic Microwave Background (CMB) maps 
with two templates constructed from the Two-Micron All-Sky Redshift Survey of galaxies. The first 
template traces the large-scale filamentary distribution characteristic of the Warm-Hot Intergalactic 
Medium (WHIM) out to ~ 90h _1 Mpc. The second traces preferentially the virialized gas in unresolved 
halos around galaxies. We find a marginal signal from the correlation of Planck data and the WHIM 
template with a signal-to-noise from 0.84 to 1.39 at the different Planck frequencies, and with a 
frequency dependence compatible with the thermal Sunyaev-Zel’dovich (tSZ) effect. When we restrict 
our analysis to the 60% of the sky outside the plane of the Galaxy and known point sources and 
galaxy clusters, the cross-correlation at zero lag is 0.064 ± 0.051/iK. The correlation extends out to 
« 6°, which at the median depth of our template corresponds to a physical length of ~ 6 — 8 /i _1 Mpc. 

On the same fraction of the sky, the cross-correlation of the CMB data with the second template is 
< 0.17 fiK (95% C.L.), providing no statistically significant evidence of a contribution from bound gas 
to the previous result. This limit translates into a physical constraint on the properties of the shock- 
heated WHIM of a log-normal model describing the weakly nonlinear density field. We find that our 
upper limit is compatible with a fraction of 45% of all baryons residing in filaments at overdensities 
~ 1 — 100 and with tem peratures i n t he ran ge 10 4 ' 5 — 10 ' 5 K, in agreement with the detection at 
redshift z ~ 0.5 of I Van Waerbeke et al.l (I2014T ). 

Subject headings: Cosmic Background Radiation. Cosmology: theory. Cosmology: observations 


1. INTRODUCTION 

The cosmological evolution of baryons within dark 
matter halos is still an open problem. While at red- 
shifts z ~ 2 the baryon fraction re siding in collapsed 
structures is close to t he cosmic mean dRauch et al.lll997t 
iWeinberg et al .1119971) . at redshifts z ~ 0 is ~ 40 — 50% 
smaller than the meas ured value (jFukugita fc Peebles I 
12004 : Shull et ah] 12012l i. Early hydrodynamical simula¬ 
tions suggested that the unidentified baryons could re¬ 
side in mildly-nonlinear filamentary structures with tem¬ 
peratures 0.01 — 1 keV and overdensities de < 10 1 2 3 the 
mean baryonic densit y, known as Warm-Hot Intergalac- 
tic Medium (WHIM) dCen fc Ostriker if 19911 iDave et, ahl 
UMI). Understanding the distribution of baryons and 
dark matter from galaxies to clusters is essential to ex¬ 
plain the formation and evolution of structure in the 
Universe. Consequently, a great observational effort is 
being devoted to detect the WHIM signature, mostly by 
searching for absorption lines due to highly ionized heavy 
elements. Some detecti ons have been reported in the 


far-ultraviolet ( 

TriDD et al.l 2008[) and in the soft X-ray 

(Nicastro et al.l 

20051: Buote et al. 2009f) bands but these 


detections are ambiguous or controversial. In some cases 
they could not be confirmed by later studies (|Yao et al.1 
l2012f ). In others it has not yet been settled whether the 
detected lines arise from the WHIM or from individual 
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galaxies or groups within the filaments ([Williams et al.1 
l2013f) . Attempts have also been made to detect signa- 
tures of the WHIM diffuse X-ray emission; IWerner et al.1 
( 2008 1 reported a detection from a filament connecting 
two clusters, with properties compatible with the WHIM. 

The baryons in this WHIM plasma do not only con¬ 
tribute to the diffuse X-ray background but they also pro¬ 
duce CMB distortions by means o f the t herm al Sunyaev- 
Zeldvich (tSZ) effect (iSunvaev fc Zel’dovich IITOTCl) . This 
effect could probe the physical state of baryons in the 
WHIM phase better than X-rays thanks to its milder 
dependence on the gas density. An alternative to find¬ 
i ng a di rect ev idence towar ds su perclu sters of galaxies 
(iMolnar fc BirkinshawllT998t IGenova-Santos et al.ll2005t) 

is to correlate CMB maps with galaxy templates tracing 
l arge scale str uc tures. This approach was first used by 
iHernandez-Monteagudo et al.l (120041 1. whose measured 
signal was dominated by clusters of gala xies. A simi¬ 
lar analysis was carried out by iSuarez-Velasquez et all 
mm an d, again, the correlation was most-likely due 
to foreground residuals and to random alignments be¬ 
tween structure in th e data a nd the matter tem plate, and 
not to WHIM. InlAtrio- Barandela fc Mucket I (120061 ) and 
I Atrio-Barandela et al.l (120081 1 we computed the power 
spectra of the thermal and kinetic SZ anisotropies 
due to t he WHIM, calcu lations t ha t were l ater re¬ 
fined by ISuarez-Velasquez, Mucket fc Atrio-Barandela - 
( 2013bli. W e used these models i i IGenova-Santos et all 
( 2009Ll2tTl3h to fit these SZ spectra to CMB data through 
a Monte-Carlo Markov Chain analysis. Although adding 
the SZ from the WHIM to the ACDM CMB power spec¬ 
trum improved the quality of the fit, this improvement 
was not statistically significant and it did not favor the 
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inclusion of a new component in the model. 

While these earlier attempts based on WMAP data 
only provided upper limits to the WHIM contribution, 
the tSZ effect has been shown to be effective at detecting 
baryonic filaments in the denser and h otter medium of 
interacting cluster pairs dPlanck Collaboration 1 [2013af) . 
This study was greatly benefited from the finer angular 
resolution and wider frequency coverage of Planck data, 
which helped to isolate be tter t h e tS Z signal from other 
components. Recently, I Van Waerbeke et al.l (j2014f) re¬ 
ported a positive correlation between mass maps recon¬ 
structed from the Canada France Hawaii Telescope Lens- 
ing Survey and a tSZ map constructed from Planck data 
with an overall significance of 6er and extending out to 
linear scales of ~ 10 Mpc. They interpreted this correla¬ 
tion as arising from WHIM at z ~ 0.4 and with tempera¬ 
tures in the range 10 5 — 10 6 K, depending on the gas bias 
and the electron density on the filaments. However, the 
measured signal could also come from gas within viri- 
alized halos, especially after the iPlanck Collaboration I 
lj2013bt) showed that the tSZ a nisotr op y exis ts down to 
masses as low as ~ 2 x 10 11 M 0 . lUrsino et al.l (l2014h sug¬ 
gested to cross-correlate X-ray and SZ maps to trace the 
WHIM distribution and found that the small-scale corre¬ 
lation would be dominated by the high-redshift WHIM, 
that again would p robe t he WHIM at early times. 

The lYan Waerbeke et al.l (I2014T) detection could indi¬ 
cate that a large fraction of the “missing” baryons re¬ 
side in a low-density warm plasma that traces the dark 
matter distribution. Equally important to this high red- 
shift observation is to probe the state of the WHIM gas 
at lower redshifts. If filaments reside in the large-scale 
dark matter potential wells, then tracers of dark mat¬ 
ter distribution at low redshifts will correlate with the 
WHIM signal both in X-ray and radio, irrespectively of 
the redshift. If the gas is confined to low mass halos 
then its tSZ will correlate with the distribution of galax¬ 
ies. In this article we compute the cross-correlation of 
the temperature fluctuations of the CMB with two tem¬ 
plates: (A) the projected matter density reconstructed 
from the distribution of gala xies within < 150/i _1 Mpc 
from the Local Group dKitaura I I2012af) to investigate 
the properties of the WHIM at redshifts z < 0.05 and 
(B) the density of the same galaxies projected along the 
line of sight to study the contribution of virialized gas 
in halos below the mass threshold of the matter tem¬ 
plate. The outline of the paper is as follows: In Sec¬ 
tion 2 we describe our data sets; we describe our proce¬ 
dure for reducing the foreground contribution of Planck 
Nomin al maps corresponding to t he March 2013 release 
(jPlanck Collaboration et al.ll2014ah and the construction 
of the galaxy and the dark matter templates. In Section 
3 we explain the filtering and the masking that we apply 
to our data. In Section 4 we present the results of the 
correlation between the templates and the CMB data; 
in Section 5 we discuss the implications on the WHIM 
physical models and, finally, in Section 6 we summarize 
our main conclusions. 

2. THE DATA 

2.1. Planck foreground cleaned maps 

As a tracer of the tSZ sig na l we u se Planck DR1 maps 
ijPlanck Collaboration et al1l2014al) . which cover a wide 


frequency range delimited by nine independent channels, 
provided by the LFI (30, 44 and 70 GHz) and the HFI 
(100, 143, 217, 353, 545 and 857 GHz) instruments. 
Apart from the cosmological signal, and the potential tSZ 
signature, these maps contain significant foreground con¬ 
tamination, mainly synchrotron and free-free emissions 
at low frequency and thermal dust at high frequency, 
as well as CO emission and zodiacal light. We care¬ 
fully clean all these contributions to obtain foreground 
cleaned maps at each frequency that will be a combi¬ 
nation of three components: CMB, tSZ and instrumen¬ 
tal noise. At the HFI frequencies, where the zodiacal 
light is more important, we use publicly-available maps 
w here this con taminant h as been rem oved, as described 
in IPlanck Collaboration et all (|2014bll . All these maps 
are produced using HealPb@ (|Gorski et al.1I2005D pix- 
elization, with resolution of N s ^ e = 1024 and 2048 for 
the LFI and HFI instruments, respectively. All Planck- 
related products used in this paper have been down¬ 
loaded from the Planck Legacy Arcliivc0. 

To perform the correction of the thermal dust emis¬ 
sion we use the maps from the Planck dust model that, 
in each sky pixel with HealPix resolution N s id e = 2048, 
contain the three parameters that define the modified 
black-body (MBB) emission law (dust-grains tempera¬ 
ture, emissivity index and optical depth) at the refer¬ 
ence frequency of 353 GHz. These maps have been ob¬ 
tained on a pixel-by-pixel basis by a y 2 -minimization fit 
of the model to the HFI an d IRAS data, as explained 
in IPlanck Collaboration et ah] (l2014d l . In order to accu¬ 
rately estimate the contribution of this emission at each 
Planck frequency we evaluate the spectral model in each 
pixel and convolve it with the passband of each detec¬ 
tor. We apply this color correction using the publicly- 
available routine hf i_colour_correction. Similarly, to 
clean the synchrotron and free-free emissions we use 
maps, delivered by the Planck Collaboration, that at 
each pixel contain the amplitude of the signal at the 
reference frequency of 30 GHz and a spectral index ob¬ 
tained by fittin g a power- law m odel to the Planck low¬ 
est frequencies (IPlanck explanatory supplement I 120131) . 
In this case a color correction using the same routine 
mentioned above, is also applied to estimate the flux 
weighted in each band. The maps of the low-frequency 
and high-frequency foregrounds are subtracted from each 
frequency map. 

The last step of the process is to clean the CO emis¬ 
sion, which is especially important in the 100, 217 and 
353 GHz channels as a result of the (1-0), (2-1) and 
(3-2) rotational transi tion li nes, respective ly. As ex¬ 
plained in the IPlanck explanatory suppl ement 1 ( 2013 ) 
and in IPlanck Collaboration et al.l ( 2014el) . the Planck 
Consortium has delivered three different types of CO cor¬ 
rection maps. The type 2 maps have the highest signal 
to noise, at the expense of having possible residuals from 
other foregrounds, but they are available only for 100 and 
217 GHz. Type 1 maps are available for the three fre¬ 
quencies, but have a much lower signal-to-noise and we 
found them not to be useful for cleaning the data. Then, 
we do not apply any correction to 353 GHz and we clean 

5 http://healpix.jpl.nasa.gov/ 

6 http: //www. sciops. esa. int/index .php?project=planck.icpage 
=Planck_Legacy_Archive 
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the 100 and 217 GHz maps using the type 2 maps, a pro¬ 
cedure that noticeably reduces the CO contamination at 
these frequencies. 

Our final foreground cleaned maps are shown in Fig¬ 
ure [U for frequencies 44 to 353 GHz. The foreground 
contamination in these maps has been clearly reduced. 
Although there are some residuals of Galactic emission 
along the plane, the signal in the rest of the sky is dom¬ 
inated by the CMB fluctuations. The power spectra 
of these maps nicely match the theoretical spectrum of 
the concordance model with the cosmological parameters 
measured by Planck. The foreground cleaning procedure 
has operated best at 353 GHz since is at this frequency 
where the Planck dust model was evaluated. The clean¬ 
ing at other frequencies is affected by possible errors in 
the MBB model or in the determination of the emissivity 
index. The residual dust contamination in the Galactic 
plane diminishes at lower frequencies. At 44 and 70 GHz 
some residual from synchrotron and free-free emission 
are visible along the plane. Although not shown here, 
the 30 GHz map was also used in our analysis. On the 
contrary, we ignored the 545 and the 857 GHz due to 
having stronger foreground residuals. 

In order to avoid sky regions affected by emission of 
diffuse foregrounds, point sources or galaxy clusters, we 
apply different masks. We first apply the Planck mask 
derived from the Comma nder-R ul er co mpo nent separa¬ 
tion methocQ (lPlanck Collaboration et aLll2014d h which 
retains 82% of the sky. Most of the masked pixels are as¬ 
sociated with point sources; a stripe along the Galactic 
plane of ~ 10° width in Galactic latitude is also removed, 
accounting for about 8% of the sky. This Galactic mask 
is sufficient to remove the residual Galactic emission that 
appear in the maps of Figure [l] To avoid regions affected 
by potentially fainter residuals we apply more conserva¬ 
tive masks. We perform the analysis using three differ¬ 
ent Planck Galactic maskfl The total sky fractions re¬ 
tained in the three resulting masks are respectively 73.3, 
55.8 and 37.1%, but we follow the practice to refer to 
these masks as the 70%, 60% and 40% masks. Finally, 
we need to eliminate possible tSZ signals from Galaxy 
clusters. Although the Commander-Ruler mask already 
removes some clusters, to be conservative we also mask 
all pixels within a radius of 0.5° around all known clus¬ 
ters with redshifts z < 0.04 and with X-ray luminosities 
Ax > 10 43 erg/s in the ROSAT [0.1-2.4]KeV band. The 
cross-correlation functions are computed on maps with 
a resolution Wide = 128. All masks are downgraded to 
this pixel resolution and, conservatively, we only retain 
pixels with values greater than 0.75. 

2.2. Density field, templates 

We will analyze two possible distributions of the miss¬ 
ing baryon fraction: (A) unbound gas following the fila¬ 
mentary structure of the large scale matter distribution, 
and (B) virialized gas in unresolved halos traced by the 
spatial distribution of galaxies. Hereafter, we will refer 
to the two templates tracing these distributions as the 
matter template M and the galaxy template M g , respec¬ 
tively. The reconstruction technique of the matter den- 

7 Mask named C0M_CompMap_Mask-rulerminimal_2048_Rl. 00 .fits 
in the Planck Legacy Archive. 

8 To be found in the HFI_Mask_GalPlane_2048_Rl. 10. f its file. 


sity field is based on a Bayesian Networks Machine Learn¬ 
ing algorithm (the Kigen-code) which self-consistently 
samples the initial density fluctuations compatible with 
the observed galaxy distribution and a structure forma¬ 
tion model g iven by second o rder Lagrangian per turba- 
tion theory (jKitaura Il2012at [Kit, aura et, al.ll2012b[ f. As 
a tracer of the matter density field we used the Two- 
Micron All-Sky Redshift Survey (2MRS, Huchra et al. 
2012). The method recovers the non-linear structures 
like knots, sheets, filaments and voids in the cosmic web 
to great accuracy. The reconstruction is performed on a 
cubic box since it requires to compute FFTs. In previous 
works we reconst ructed the density fie ld using boxes of 
side 160fe~|Mpc dKitaura et all [2(112 Hi and 180/i _1 Mpc 
(|Hefi et all 1201311 and found compatible results, proving 
that the boundary effects are unimportant. In this work, 
the reconstructed density field on the sphere includes 
modes with wavelengths up to 110 — 130/i _1 Mpc. For 
our theoretical estimates we take the average depth of 
the reconstructed density field to be 90/i _1 Mpc. The 
template built through this reconstruction technique is 
produced using a HealPix resolution Wide = 128. The 
matter distribution shows well defined filaments, charac¬ 
teristic of the mildly non-linear regime that we assume 
to be a good tracer of the WHIM. As the reconstruction 
also traces the galaxy distribution, it is sensitive to the 
hot halos surrounding galaxies that have been recent ly 
detected (|Anderson fc Bretnnan |[20TH : I Dai et al.ll2012lf . 

Low-mass galaxy clusters or galaxy groups below the 
threshold mass of our previous tem plate co uld also pro¬ 
du ce a detectable tSZ anisotropy (lPlanck Collaboration I 
l2013bll . The pixels at Wide = 128 have an angular size 
of ~ 0.5° x 0.5° and if the virialized and unbound gas 
were to follow the same spatial distribution we would 
not have enough angular resolution to disentangle a po¬ 
tential signal associated to collapsed low-mass systems 
from the signal of the unbound WHIM gas. Most likely 
the virialized gas detected by Planck is better traced by 
the galaxies themselves. Under this assumption, we can 
check the amplitude of this contribution using a tem¬ 
plate of the projected density of galaxies. To make the 
comparison more meaningful, we project along the line 
of sight the same sample of 2MRS galaxies that went 
into the reconstruction of the matter density field de¬ 
scribed above. By construction, this galaxy template 
traces the hot circumgalactic halo. At 10/i -1 Mpc a pixel 
of half degree side corresponds to 80/i _1 Kpc, slightly 
larger than the halo but at 60/i _1 Mpc, the average 
depth of our template, it corresponds to ~ 500/i _1 Kpc, 
much larger than the vi rial radius of the detected halos 
ij Anderson fe Bregman I [201 111 . Therefore, this template 
could potentially also trace the gas outside the circum¬ 
galactic halo. 

3. DATA PROCESSING 

The tSZ signal due to WHIM has an amplitude of 
about few pK but the expected cross-correlation with 
a matter template out to ~ 9 0 — 1 20fe~ 1 M pc is ex¬ 
pected to be ~ 0.03 — 0.3 pK (jSuarez-Velasquez et all 
!2013aT) . At any frequency u, the error bar on the cross¬ 
correlation of the CMB data T{y) with a matter tem¬ 
plate M, (T(u)M)(6), is dominated by the intrinsic CMB 
anisotropy. In order to obtain a statistically significant 
detection, the data needs to be processed further. There 
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are two possibilities: (1) use the 217 GHz channel, the 
frequency of the tSZ null, to subtract off the intrinsic 
CMB signal while preserving the tSZ anisotropy and (2) 
construct a Compton parameter y-map combining all fre¬ 
quencies to remove all compo nents from th e data except 
the tSZ, as it was done in iPlanck Collaboration et ahl 
(|2014sD . Subtracting the 217 GHz channel from the other 
maps will increase the noise and add the 217 GHz fore¬ 
ground residuals to those remaining at other frequen¬ 
cies but, when constructing a y-map the frequency infor¬ 
mation of the tSZ effect will be lost. This information 
could be crucial to distinguish correlations due to ran¬ 
dom alignments or stripes present in the data from the 
true effect due to the WHIM. Therefore, in this article we 
will use method (1), i.e. we will subtract the correlation 
function at 217 GHz from those at other frequencies. 

In Figure [2] we plot the power spectra of the subtracted 
maps. We only analyzed maps with HealPix resolution 
fVsMe = 128. Effectively, all the maps have been down¬ 
graded to the same resolution by the pixel window, which 
dominates over that of the beam. When subtracting the 
217 GHz from the maps at other frequencies, we are effec¬ 
tively removing the intrinsic CMB signal. At multipoles 
t < 20, the power spectra scale as Ce ~ £~ 2 , compat¬ 
ible with being dominated by marginal CMB residuals, 
instrumental (1//) noise and foreground residuals, while 
at l > 60 the spectra are compatible with white noise. 
Maps could be correlated even if there is not any WHIM 
signal present. To reduce these spurious correlations, 
maps can be filtered by removing the lowest multipoles 
up to any given t cnt . Since at N s id e = 128 the highest 
multipole is £ max ~ 256, removing modes up to £ cu t = 60 
could erase a significant part of the WHIM signal. As 
a compromise, we set £ cut = 20. Finally, to reduce the 
transfer of power from galactic residuals to multipoles 
at all scales, before computing the multipole expansion 
of the original map, 20% of the sky around the galactic 
plane is masked. 

In Figure [3] we show the 353-217 GHz difference CMB 
map (top row), the matter template (middle) and the 
galaxy template (bottom) before and after the lowest 
multipoles have been removed. If initially the CMB data 
is dominated by noise and foreground residuals, cutting 
the lowest multipole reduces both contributions. Stripes 
due to Planck scanning strategy are still present and, 
with the smaller noise levels, they are more clearly seen 
in the filtered map. Filtering introduces aliasing between 
different multipoles, giving rise to significant tempera¬ 
ture anisotropies next to the galactic plane; in addition, 
fringing is also present. Nevertheless, these artifacts do 
not pose a serious problem. In fact, we can compute the 
CMB-template cross-correlation with different masks, re¬ 
moving from 20% to 60% of the sky, to test the effect of 
these systematics. 

For consistency we also remove the lowest multipoles 
of the templates. The matter template is shown in the 
middle row of Figure [3] before (left panel) and after (right 
panel) filtering. Since the template has little power at 
large scales (see Figure 2 of Suarez-Velasquez 2013a), 
the filtered template has very similar structure to the 
original template and the filaments show a very similar 
distribution. The template of the projected density of 
galaxies before (left) and after (right) filtering is shown 
in the bottom row of Figure [3] This template has a sim¬ 


ilar large structure as the matter distribution but is less 
smooth at small scales as one could expect of the spa¬ 
tial distribution of clustered halos. Again, notice that 
filtering the lowest multipoles does not affect the overall 
distribution of galaxies. 

A visual comparison of the filtered templates of mat¬ 
ter (right middle panel of Figure^]) and projected density 
of galaxies (right bottom panel) shows that they have a 
very similar large scale structure, as one could expect 
since the galaxy template contains all galaxies that were 
used to reconstruct the matter distribution. Their main 
difference is at small scales where the matter template 
is a smoother function. Then, one can assume that the 
matter template would trace better the unbound WHIM 
gas residing in the large scale dark matter potential wells 
while the galaxy template would trace the circumgalactic 
gas in virialized halos. However, as discussed in Sec. 2.2, 
both templates could trace both gas components, proba¬ 
bly with different relative weights, since the matter tem¬ 
plate includes the contribution from the hot halos around 
galaxies and for most galaxies the angular scale sub¬ 
tended by the virial radius of those circumgalactic ha¬ 
los is smaller than the pixel size at Healpix N s i(j e = 128 
resolution. 

To summarize, our pipeline performs the following 
steps on the foreground cleaned Planck Nominal maps: 
first, the galactic plane, point sources and clusters are 
masked out. Second, CMB data and templates are 
transformed to spherical harmonic space and all mul¬ 
tipoles with £ < 20 are removed. Next, all maps are 
downgraded to a HealPix resolution of Wide = 128, the 
resolution of the matter template. Templates are nor¬ 
malized to zero mean and unit variance so the cross¬ 
correlation is given in temperature units. Finally, the 
cross-correlations C v = {T v * M) are calculated at each 
frequency v using three different masks that preserve 70, 
60 and 40 per cent of the sky. 

4. RESULTS 

In order to remove the contamination from the primor¬ 
dial CMB, we subtract the 217 GHz correlation function 
from those calculated at all other frequencies. If the 
cross-correlation (C v — G 217 ) is dominated by the tSZ 
effect due to the WHIM traced by the template, then 
it must scale as (C v — G 217 ) ~ G(v), reflecting the fre¬ 
quency dependence of the tSZ effect G(is). In this case, 
the reduced cross-correlation ( TM) = ( C v — C 217 )/G{v) 
between a template and the CMB data must be in¬ 
dependent of frequency and proportional to the cross¬ 
correlation of the template with the Comptonization pa¬ 
rameter Y c , i.e., {TM) oc (Y C M). This reduced correla¬ 
tion function is given in Figure[4j the top row corresponds 
to the correlation with the matter template (Figure [3] 
middle right) and the projected galaxy density template 
(Figure [3] bottom right). Panels (a), (d) correspond to 
the correlation with LFI frequency maps and (b), (e) to 
the correlation with HFI frequency maps computed on 
60% of the sky; (c) and (f) is the combination of all fre¬ 
quencies for different fractions of the sky. The tSZ is 
negative in the Raleigh-Jeans part of the spectrum and 
is positive in the Wien region. The frequency depen¬ 
dence of the tSZ effect varies from G(30 GHz)~ —2 to 
G(143 GHz)~ -1 and G(353 GHz)~ 2.2. Then, if the 
cross-correlation of the templates and Planck foreground 
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cleaned data corresponds to a tSZ contribution, it must 
be negative for the frequencies below 217 GHz and pos¬ 
itive above it. This is what it actually happens. Panels 
(a), (b), (d) and (e) of Figure 0] show that all correlations 
have the same sign once divided by G{v ), indicating that 
the they are consistent with being due to tSZ contribu¬ 
tion since their amplitude is similar in all channels and 
have the correct sign and, although none of our measure¬ 
ments is statistically significant, the correlations with the 
matter template in the top row are a factor three larger 
than those with the projected galaxy density in the bot¬ 
tom row. At the origin, the signal-to-noise ratio (SNR) 
is of order unity, being somewhat larger at angular scales 
(9~2 —5°. 

In Figure [4] error bars were computed from the rms 
dispersion of the correlation of 100 random rotated tem¬ 
plates with the CMB maps. To explore the sky homo¬ 
geneously, we select a random set of N s j m pixels drawn 
from a random distribution of all the pixels in the map 
and compute their angular coordinates (l,b). Then, we 
compute the Euler angles that align the north galactic 
pole with the direction (l, b) of the randomly chosen N s i m 
pixels and rotate the template. Since our template re¬ 
constructs the full sky, the area on which the correlation 
is computed is always defined by the mask, and there¬ 
fore is the same in all the rotations. By using the same 
data to compute the correlation and its error bars we 
include all contributions due to foreground residuals and 
1 // instrumental noise. 

In panels (c) and (f) of Figure 01 we plot the correla¬ 
tion functions after combining the six LFI+HFI channels. 
For comparison, we plot the results for the three masks 
used in our study. To combine different frequencies, for 
each angle we make a weighted average of the correlation 
functions for different frequencies: 


<TM) av (0) 


J2,{TM)(d,y l )w(e i Vi) 

E * w ( 6 i v i) 


(1) 


where the sum runs over all the frequencies considered, 
and w(6,i , i) = 1/<t( 0, zzj) , with a(9,i/i) the rms disper¬ 
sion calculated at position 9 in all the rotations of map z/,. 
We calculate the error bar associated to this estimator by 
accounting for possible correlations between frequencies: 

cr((TM) av (9)) = 1 = , 

\JHi Ej w{6, Vi)(TM)ij(9)w(9, uj) 

( 2 ) 

where ( TM)ij(9 ) = (( TM)(9,Vi)(TM)(9 1 Uj )) represents 
the covariance between frequencies i/j and vj at position 
9. We found that the results for difference frequencies 
are strongly correlated and for this reason the statistical 
significance does not improve much by averaging differ¬ 
ent frequencies. In particular, at zero lag and for the 
60% mask, where the covariance is typically between 0.75 
and 0.85, we obtain a SNR= 1.26, while the signals for 
individual frequencies have SNR between 0.84 and 1.39. 
Our overall statistical significance is never larger than 
SNR~ 2. After averaging over all frequencies, the corre¬ 
lation with the matter template, measured on 60% of the 
sky, was (TM)( 0) = 0.064 ±0.051 fiK. Taking this result 
as an upper limit we obtain (TM)( 0) < 0.11 fxK at 68% 
C.L. and < 0.17 fiK at 95% C.L. Even if the result is not 
statistically significant, we find that the cross-correlation 


is zero at 9 ~ 6 — 8 degrees. The average depth of our 
template is ~ 60/i _1 Mpc. At this depth, the zero cross¬ 
ing of the correlation function corresponds to a linear 
scale of ~ 6 — 8 /i _1 Mpc. This would be the averaged 
projected linear size of the WHIM filaments traced by 
our matter template. Thi s size is reassur ingly consi stent 
with the scale found by IVan Waerbeke et all (l2014h and 
could indicate that the effect is real but our template is 
not deep enough to produce a significant detection. 

Comparison of the cross-correlation using different 
galactic masks allow us to test the effect of fringing and 
power leakage of the galactic foregrounds into the galac¬ 
tic poles due to our filtering scheme. The oscillation pat¬ 
tern at large distances of the correlation functions shown 
in Figure [4] is partly produced by this residual fringing. 
In addition, the number of pairs varies with a scale of 
~ 6° — 10°, the mean separation between filaments and 
also the typical angular size of filaments, increasing and 
decreasing the sample variance component of the error 
bar accordingly. Notice that when we restrict the analy¬ 
sis to the cleanest regions of the sky, located around the 
Galactic Poles, the statistical significance of our results 
increases because even if the error bars are larger due 
to a bigger sampling variance the signal increases even 
more. 

If the galaxy template traces preferentially the gas 
stored in halos, it can be used to constrain the contri¬ 
bution due to the halos below the mass threshold of the 
galaxy clusters excised from our analysis. The cross¬ 
correlation between the galaxy density template and the 
foreground cleaned maps at zero lag in 60% of the sky is 
also compatible with zero: (TM g )( 0) = 0.022±0.028 /iK. 
We obtain an upper limit of (TM g )(0) < 0.078 at 
the 95% confidence level. As we have argued, the mat¬ 
ter and the galaxy template trace both the bound and 
unbound gas components, probably weighing differently 
each contribution. Comparing with the previous result, 
the galaxy template provides a more restrictive upper 
limit that could reflect that the template does not trace 
the WHIM fully. 

5. CONSTRAINTS ON WHIM MODEL PARAMETERS 

We can use the 95% upper limits given 
in the previous section to constrain the 
phy sical state o f the gas in fila ment s. In 
iSuarez-Velasquez. Miicket fc Atrio-Barandelal (l2013bf) 
we described the distribution of baryons in the WHIM 
as a network of filaments following a log-normal dis¬ 
tribution function. In our model baryons follow the 
dark matter distribution except at small scales, where 
baryon perturbations are damped by shock heating 
(|Klar fc MiicketlfeOlOf) . At any given redshift, the cut-off 
scale L cnt is determined by the condition that the linear 
velocity perturbation v(x, z) is equal to or larger than 
the sound speed c s = (/cb7igm(z)/top) 1//2 , with m p the 
proton mass and k b the Boltzmann constant. Then, the 
cut-off scale L cut can be parametrized in terms of the 
mean Intergalactic Medium temper ature, Ttgm. While 
the latter changes little with redshift (jTittlev fc Meiksinj 
12001 . L cut « L^l + z) 1 / 2 . At Tigm = 10 4 K, the cut-off 
scale today would be Lq = 1.7 h~ l Mpc (for details see 
Suarez-Velasquez, Miicket & Atrio-Barandela, 2013b). 

To compute the tSZ effect generated by the free elec- 
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trons in the WHIM, we need to specify the equation 
of state. For shock -he ated gas, the phase diagrams ob¬ 
tained by iKang et all (|2005j ) from numerical simulations 
can be fitted by log 10 (T e /10 8 .K) = —2[log 10 (3.5 + x 6 )] -1 
with ft = a + x^ 1 and x = n e /fiB the electron den¬ 
sity in units of the mean baryon density. In addition 
to the two parameters Tigm and a, WHIM temperature 
anisotropies depend on erg, the amplitude of the matter 
power spectrum on a sphere of 8 ft -1 Mpc, and on the 
fraction of free electrons in the filaments. In Figure [5] 
we present the radiation power spectra of the WHIM 
temperature anisotropies for different model parameters; 
we assume that all the baryons in the WHIM are ion¬ 
ized. In these mo dels we fixed the Pla nck measured value 
of erg = 0.83 dPlanck Collaboration 1(20 14f) . This corre¬ 
sponds to a WHIM baryon fraction of 43% to 48% for 
Tigm = 10 3 - 6 , 10 4 K, respectively. The spectra in panel 
(a) of Figure [5] represent the full contribution, obtained 
by integrating up to redshift z up = 3. The contribution of 
filaments from earlier times is negligible. Spectra are pre¬ 
sented in pairs, corresponding to Tigm = 10 3 6 K (solid) 
and 10 4 K (dashed lines), respectively. From top to bot¬ 
tom, curves correspond to a equation of state parameter 
a = 2.0, 0.9, respectively. In panel (b) of Figure [5] the 
power spectra follow the same convention. In this case 
the integration is restricted to z up < 0.03, similar to the 
scale probed by our template. 

Our formalism allows us to predict the cross¬ 
correlation of the spatial distribution of the WHIM with 
the tSZ distortion generated by it. As indicated in sec¬ 
tion the matter template describes the matter dis¬ 
tribution in the mildly non-linear regime. If the ionized 
gas follows the matter and its equation of state is poly¬ 
tropic, T e oc ni 7 1 \ then the pressure profile will be 
P e oc n/, where 7 = 1 corresponds to isothermal gas, 
characteristic of the shocked regions and 7 = 5/3 to adi¬ 
abatic monoatomic gas. From the template M oc n e we 
have constructed pressure templates with 7 = (1,5/3). 
We verified that the cross-correlation at zero lag var¬ 
ied by less than 20 % for adiabatic indices in the above 
range, indicating that if the template describes the dis¬ 
tribution of the ionized WHIM, then it will also trace 
its tSZ contribution. Then, if a matter template, nor¬ 
malized to zero mean and unit variance, traces the tSZ, 
the cross-correlation would be (T M) ( 0 ) = ctwhim where 
°whim = + l)C/ VHIM / 47 r; therefore, the results 

presented in Figure [5] give the expected correlation for 
each set of model parameters. Additionally, we need to 
consider that the actual realization of the WHIM in the 
sky is not exactly the mean power spectrum since our 
template is just one realization of the sky. This gives rise 
to an uncertainty known as cosmic variance and given by 
ecv = J2 \/2(21? + l)Ce/4n. 

In Figure [B] we present the amplitude of the cross¬ 
correlation at zero-lag, computed using the upper limit 
of integration z up = 0.03, similar to the depth of our 
template with the upper limits obtained from the data. 
In panels (a) and (b) of Figure [ 6 ] the model corresponds 
to Tigm = 10 36 ,10 4 K, respectively. Plots follow the 
same convention; the thick solid blue lines represent the 
theoretical prediction and the surrounding shaded green 
area the ler cosmic variance uncertainty. The horizon¬ 


tal (red) lines correspond to the 95% confidence limit 
given by the measured cross-correlation with the matter 
and the galaxy templates. The dashed (violet) line rep¬ 
resents the correlation for an upper limit of integration 
z U p = 0.05 to show that by increasing the depth of our 
template we could probe the model very effectively. Since 
our error bars are dominated by the random alignments 
of our filaments with large scale structures in the CMB 
data, they will not increase significantly if our template 
was to trace the matter distribution out to, for example, 
that redshift but the signal will rise by almost a 50%. If 
we were not to have a detection then the constraints on 
model parameters would be much stronger. 

Using a grid of models in the range as = 
[0.7,1], a = [1,4] and log 10 (Ti GM ) = [3.4,4.2] 
we have verified that the power spectra scales with 
cosmological and model parameters as ctwhim oc 
(logTiGM) _1 ' 4±0 ' 2 a 1 ' 7±02 c r 8’ 3±0 ' 1 UB. Since as and Hb 
are fixed by observations the results of Figure|B]show that 
the equation of state parameter a is the most important 
parameter of the model. In the observationally allowed 
range, the variation of the model prediction with Tigm is 
rather weak. Including cosmic variance and considering 
the upper limit at the 95% confidence level of the matter 
template, a < 1.6 and < 1.9 for Tigm = 10 36 K and 
10 4 K, respectively. The bound is much more stringent 
when we impose the upper limit derived from the galaxy 
template: (TM g }( 0) < 0.08 /rK: a < 1.1 and a < 1.3, 
respectively for the two values of Tigm considered. 

The upper bound in a can be translated into a con¬ 
straint on the temperature of the gas at different den¬ 
sities using the equation of state derived from the re¬ 
sults of Kang et al (2005). Figure [7] defines the range 
of the temperature of the shocked-heated WHIM as a 
function of its density contrast. The solid line corre¬ 
sponds to a = 1.0 and the dashed lines correspond to 
the 95% confidence level upper limit from the matter 
template, a = 1.9 (upper line) and the same limit from 
the galaxy template, a = 1.3 (lower line). The shaded 
green area shows the region excluded by the more strict 
upper limit from the galaxy template that is still allowed 
by the matter template. For comparison we also plot 
th e same temperature fun ction for the numerical results 
of lCen fe Ostrikerl (|2006f h For instance, at overdensity 
x = 1 the temperature of the electron gas is T e ~ 10 4 8 K, 
while for x = 100 it is T e < 10 7 5 K. These tempera¬ 
ture ranges are in agre e ment with the results of simu - 
lations (see IKang et all (|2005l k iCen fo Ost.riker I (|2006lh 
IBertone et aD (IMl l. Still, at Tigm = 10 36 - 10 4 K the 
fraction of baryons in the WHIM medium varies from 
43% to 48%, providing enough room to accommodate 
the missing baryon fraction in the local Universe. 

6. CONCLUSIONS 

We have searched for the tSZ signature of the WHIM in 
Planck Nominal maps. To that purpose we have cleaned 
the frequency maps of Planck first data release from fore¬ 
ground emissions. These maps were correlated with a 
template reconstructed from the observed galaxy dis¬ 
tribution in the 2MRS that describes the matter den¬ 
sity field out to « 90ft -1 Mpc that traces the large- 
scale filamentary distribution of the WHIM and with a 
galaxy template that traces the hot circumgalactic ha¬ 
los. Our results showed no evidence of a WHIM-induced 
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tSZ temperature anisotropy to be present at any signif¬ 
icant level or of any detectable contribution from viri- 
alized gas within low-mass halos. We set upper limits 
of < 0.08 fj,K and < 0.17 /iK at 95% C.L. using the 
results from the galaxy and matter templates, respec¬ 
tively. We translated these upper limits into constraints 
on the p roperties of the WHIM using the m od el de- 
scribed in lSuarez-Velasquez. Mticket fe Atrio-Barandela I 
(|2013bD . We found that the temperature of the shock- 
heated WHIM gas was in the range T = 10 4 ' 5 — 10 7,5 K at 
overdensities in the range 1 — 100. The cross-correlation 
of the matter template and the Planck CMB data goes to 
zero at 9 ~ 6°. At the median depth of our template, this 
angular scale would be compatible with filaments having 
an averaged projected size of 6 — 8/i _1 Mpc. The size of 
filaments and the temperature of their gas is in excellent 
agreement with the scale and temperature range found 
by Van Waerbeke et al (2014), indicating that the matter 
template could be tracing the WHIM distribution even 
though it is not deep enough to produce a statistically 
significant signal. 

The search for the missing baryons is not restricted 
to the CMB. Evidence of the distribution of this mat¬ 
ter component has been obtained via UV-absorption line 
studies with FUSE and HST-COS (Shull et al. 2012) 
but these observations probe the coldest fraction of these 
baryons. Our results imply that most baryons would 
have temperatures T > 10 4 5 K, in agreement with hy- 
drodynamical simulations. This phase could be detected 
through highly ionized C, N, O, Ne, and Fe ions. In this 
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respect, the proposed Athena X-ray mission (Kaastra 
et al 2013) could provide very useful information about 
the WHIM parameters. With WHIM gas temperatures 
in the range T = 10 4 5 — 10 7 5 K, the filaments in our 
model can accommodate up to 48% of all the baryons. 
Since the tSZ effect is proportional to the electron pres¬ 
sure, if X-ray observations were to find that the average 
WHIM is larger than the limits given in Figure [7| then, 
to make compatible the tSZ upper limit with a higher 
WHIM temperature, the baryon fraction stored in fila¬ 
ments must decrease in a proportion similar to the tem¬ 
perature increment. Then, by combining observations 
of the WHIM at different wavelengths and different red- 
shifts, we will measure better the physical properties of 
the WHIM, provide better constraints on the model pa¬ 
rameters, improve our understanding of the physical pro¬ 
cesses undergone by the baryons during the formation of 
galaxies and the emergence of large scale structure, and 
could separate the contribution of gas in unresolved halos 
from that of the WHIM. 
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Fig. 1.— Foreground cleaned Planck Nominal maps at different frequencies. From top to bottom and left to right, frequencies 
are 44, 70, 100, 143, 217 and 353 GHz. A detailed explanation of the cleaning methodology is given in section [2d] Apart 
from some small foreground residuals along the Galactic plane (which are positive in the case of the thermal dust emission and 
negative for the synchrotron emission), these maps are dominated by primordial CMB. 
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Fig. 2. Power spectra of the foreground cleaned Planck Nominal maps, after dividing by the window function of each frequency 
and subtracting the 217 GHz channel. The increase in power at low multipoles could be produced by a combination of detector 
1 // correlated noise and foreground residuals. To remove this effect we filter the cleaned maps by removing all multipoles with 
t < 20, as explained in section [3] 
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Fig. 3.— Comparison of the original CMB data and templates with their filtered counterparts obtained by removing all 
multipoles with l < 20. Top row: original Planck Nominal map at 353 GHz minus 217 GHz (left) and filtered map (right). The 
mask applied to the data removes 20% of the sky. Middle row: original (left) and filtered (right) matter template. Bottom row: 
original template of projected density of galaxies (left) and filtered (right) template. 
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Fig. 4.— Cross-correlation of filtered matter (top row) and galaxy density (bottom row) templates with the foreground cleaned 
Planck Nominal maps divided by the frequency dependence of the tSZ effect; at the different frequencies the cosmological signal 
was removed by subtracting the 217 GHz channel. Solid lines represent the data and dashed lines the error bars. In (a,b) and 
(d,e) we plot the correlations with the LFI and the HFI channels, respectively, calculated using the mask that retains 60% of 
the sky. In (c,f) we plot the average correlation of all frequencies obtained after masking different fractions of the sky. The 
percentage indicates the fraction of the sky used. 
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Fig. 5.— Power spectra of the tSZ due to the WHIM for different model parameters. As indicated the spectra in panel (a) 
includes the contribution at all redshifts, while the spectra in panel (b) only includes the contribution up to z up = 0.03. Solid 
and dashed lines correspond to Tigm = 10 3 ' 6 K and 10 4 K, respectively; neighboring graphs, given in the same color, correspond 
to the same value of the equation of state parameter a, in decreasing amplitude from top to bottom. 
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Fig. 6.— Comparison of the measured and the predicted amplitude of the cross-correlation of the matter template that traces 
the matter distribution out to a up = 0.03 and the tSZ anisotropy due to the WHIM. The horizontal red lines give the upper 
limit of the measured correlation at zero lag at the 95% confidence levels derived from the matter template and the galaxy 
templates. The blue solid line and the shaded green area represent the expected cross-correlation and the 1 <t cosmic variance 
uncertainty for different model parameters, considering a maximum redshift z up = 0.03 for the matter density field traced by 
the template. The dashed violet line represents the predicted cross-correlation for z up = 0.05. In panel (a) and (b) the IGM 
temperature is Tigm = 10 3,6 K and 10 4 K, respectively. 
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x = n e /n B 


Fig. 7.— The blue solid line represents the temperature of the shock-heated WHIM for different overdensities in the Kang et 
al (2005) numerical fit for an equation of state parameter a = 1. The blue dashed lines correspond to a = 1.3 (bottom) and 
a = 1.9 (top), the 95% upper limits from the galaxy and matter templates, respectively. The shaded green area shows the 
region excluded by the more strict constraint of the galaxy template that is allowed b y the matter template. For comparison, 
the dot-dashed red line represents the temperature-overdensity relation from ICen fc Ostriker I (1200611 . 


15 












